Impurity effects in the quantum Kagome system ZnCu 3 (OH) 6 Cl2 ? 
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Motivated by the recent experiments on the new spin half Kagome compound ZnCu3(OH)6Ci2, 
we study a phenomenological model of a frustrated quantum magnet. The model has a spin liquid 
groundstate and is constructed so as to mimic the macroscopically large quasi-degeneracies expected 
in the low- lying energy structure of a Kagome system. We use numerical studies of finite size 
systems to investigate the static as well as the dynamical response at finite temperatures. The 
results obtained using our simple model are compatible with a large number of recent experiments 
including neutron scattering data. Our study suggests that many of the anomalous features observed 
in experiments have a natural interpretation in terms of a spin-j defects (impurities) coupled to an 
underlying Kagome- type spin liquid. 

PACS numbers: 75.50.Ee, 75.10.Jm, 75.40.Gb 
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The quest for spin liquids in real materials has been 
in the forefront of condensed matter physics since An- 
derson's suggestion of resonating valence bond states as 
a possible ground state for the antiferromagnetic insu- 
lating phase seen in high Tc superconductors. Typi- 
cally one expects spin liquid ground states in geometri- 
cally frustrated systems with a high degree of frustration 
like in the Heisenberg antiferromagnet on pyrochlore lat- 
tices, checkerboard lattices and the Kagome lattice 
The Heisenberg antiferromagnet on the Kagome lattice 
is particularly interesting in that numerical calculations 
show the existence of a special kind of spin liquid ground 
state with a macroscopic quasi degeneracy [2j. However, 
not much is known about the dynamical properties of 
these spin liquids. Interest in the physics of Kagome 
lattices has been boosted by the synthesis of a novel 
paratacamite compound ZnCu3(OH)gCi2, which is con- 
sidered as a faithful realization of a spin-i Kagome sys- 
tem Q . A spate of recent experiments 0, IE Hi 0] have 
however, generated a huge debate about whether the ob- 
servations are actually due to intrinsic Kagome physics, 
or whether other aspects like impurities or Dzyaloshinski- 
Moriya interactions 8] may need to be taken into account. 

Motivated by the current lack of perfect stoichio- 
metric control in the synthesis of ZnCu3(OH)gCl2, and 
its impurity-like low temperature behavior of the static 
magnetic susceptibility, we introduce a schematic model 
which toys with the idea that the anomalous features seen 
in experiments may arise from spinful impurities dressed 
by an underlying spin liquid with an energy spectrum of 
the Kagome type. Our mean-field model explicitly mim- 
ics the principal feature of the Kagome antiferromagnet: 
a high degree of frustration leading to a macroscopic de- 
generacy of the ground state. Due to its mean-field na- 
ture, our model is much easier to study than the full 



Kagome system. We begin by considering the following 
hamiltonian, 



H = 
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where, S represent the spin half operators and Jy the 
antiferromagnetic interactions between two spins. Since 
all the spins interact with each other, the model is max- 
imally frustrated. For Jjy = J for all the model 
is exactly soluble since the hamiltonian can be written 
as H = J/2[Sf ot — 3N/4] and all states with a given 
value of total spin S tot are degenerate. For N even, the 
ground state is a bunch of degenerate singlets with a 
gap to the degenerate triplets given by A = J. For N 
odd, the ground state is described by degenerate spin 
\ states with a gap to the quadruplet spin | given by 
A = 1.5 J. In order to mimic the quasi degenerate na- 
ture of the low energy sector of the Kagome spin sys- 
tem, we add a small disorder term to the magnetic in- 
teractions Jij = J(l + 6ij) where Sij << 1 are nor- 
mally distributed random numbers. This small change 
in couplings is sufficient to lift the exact degeneracy of 
the states, but not large enough to destroy the energy 
separation between the low-lying sectors. The number 
of quasi degenerate states of the ground-energy sector is 
macroscopic and given by A even w 1.68"" for N even, 
and A dd ~ 1.74" for N odd. This resembles the scaling 
seen in the Kagome lattice where A even « 1.14" for N 
even and A ot id ~ 1.16" for N odd. Thus, in the N even 
(odd) case both models feature a quasi-degenerate and 
macroscopically large lowest energy sector with S = 
(S = 1/2), and a first excited sector of similar character- 
istics with S — 1 (S = 3/2). We solve the model using 
exact diagonalization methods 0, U3] and calculate the 
static spin susceptibility and the local dynamic suscepti- 
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FIG. 1: (Color Online) Typical behavior of the static suscep- 
tibility as a function of temperature T for N odd (full line) 
and N even (dashed line). The inset shows the typical tem- 
perature variation of the normalized magnetization in a field 
for N odd. The finite intercept and the convex shape of the 
curve are in qualitative agreement with the experiments of 
Refill 

bility in the presence and absence of an external magnetic 
field. These calculations are done at finite temperatures 
which restricts the size of the systems that we can access 
numerically. Nonetheless, as we will show below, finite 
size effects are rather minimal. This is partly due to the 
mean field nature of the model, and, also to the fact that 
the computed observables are obtained as direct averages 
over realizations of disorder [jj [nj] • 

Regarding the current debate on whether the observed 
experimental behavior is intrinsic to the putative spin liq- 
uid state of a Kagome lattice or due to impurity effects, 
we shall show below that many experimental observations 
can be naturally explained as due to strongly interacting 
many-body spin-i states that are present only in N odd 
systems. We shall argue that this leads to the interpreta- 
tion that low energy features seen in current experiments 
on the paracatamite are controlled by impurities. 

To begin with the discussion of our toy model re- 
sults, we first consider the uniform static susceptibility 
per spin, x(T). This quantity is easily computed from 
the magnetic moment m(h) for small fields i.e., x—mjh^ 
where m is the magnetic moment per spin at a temper- 
ature T and h is a small external applied field. Due to 
the mean-field nature and the translationally invariant 
spin liquid ground state of our model, the local and the 
global magnetic moment are the same. Consequently, the 
ensuing analysis holds for local as well as global suscep- 
tibilities Q. From Fig. [TJ we see that x(T) is strongly 
dependent on whether the total number of spins is even 
or odd. This can be easily understood from the fact 
that even systems have a singlet ground state (with other 
singlets nearly degenerate with the ground state) and a 
well defined gap to the first triplet. This gap leads to 
X(T) (X exp(-A/T) -> as T -> 0. On the other hand, 
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FIG. 2: (Color Online) The finite size scaling for the sus- 
ceptibility show the saturation of \ as T — > in a field. The 
dotted green line represents \ = 1/4T. The inset shows the 
variation of the coefficient a defined in Eq.(2) with system 
size. 

systems with TV odd have a multitude of low-lying spin 
doublet states that are quasi-degenerate with the ground- 
state. Thus, as T — > 0, this results in a Curie like behav- 
ior of the susceptibility 

«m»f (2) 

The inset of Figj2] shows the scaling of the Curie coeffi- 
cient a (TV). Interestingly, the finite size scaling reveals 
that, despite the macroscopic number of low- lying spin 
doublet states, a(N) = l/N. Thus, the strength of the 
divergent part of x(T) corresponds exactly to the contri- 
bution of a single effective spin-^, independently of N. 
In other words, since x(T) is the susceptibility per spin, 
our result implies that for systems of any size N (odd) 
the strength of the Curie- Weiss contribution of the total 
spin susceptibility N\(T) is exactly 1. This is a clear 
indication of a spin-i many-body state emerging in odd- 
N systems, which can be also thought of as a missing 
spin-i (defect) in an otherwise S — spin liquid. The 
existence of this effective " dangling" spin is further clari- 
fied by the behavior of x(T) at small finite fields h. From 
FigOl we see that the Curie divergency observed in odd 
N systems is cut-off by h, and x becomes temperature 
independent for T < h. This is in agreement with the re- 
sults reported in RefH, where a field of 0.2T was applied 
and the saturation of x(T) below 0.2K was observed. At 
intermediate and high temperatures, x(T) is almost the 
same for systems with odd/even N modulo finite size 
scaling effects. At high enough temperatures T > J, we 
recover the usual Curie- Weiss behavior of the suscepti- 
bility x(T) °c {T + lew) -1 - As expected for a mean 
field model where every spin has N-l neighbors, we find 
that the Curie- Weiss temperature which fixes the scale 
of the antiferromagnetic exchange in the system is given 
by Tew ~-0.3(N-l) J where the factor 0.3 is due to quan- 
tum fluctuations and the inherent frustration. To com- 
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pare the scale of the Curie- Weiss temperature predicted 
by the mean field model with that of real experiments, we 
can use the experimental value of J=170K. For Kagome 
lattices, since the number of nearest neighbors is 4, we 
obtain |X"cw| =1.2J (~204K) which is of the same order 



as the experimental estimate of around 270^320K. 

We now discuss the behavior of dynamical quantities. 
We use the following spectral decomposition to calculate 
the imaginary part of the dynamic susceptibility Imx = 
X": 



x"H = 



M N 



ZNM 



2^2. \^ m \S> m )V^ - E™ + £^)exp(-^)[l - exp(-/3c)] 



m=l i—1 vy. 
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e — o 7 sites 
h — a 9 sites 
a— * 1 1 sites 
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FIG. 3: (Color Online) Fhe T=0 dynamic susceptibility x"(w) 
as a function of frequency id in units of J for systems with N 
odd and N even. 



where Z is the partition function and M the number of 
disorder realizations that are used to average x- I n Fig- 
[3] we plot x" f° r both even and odd systems. Note that, 
as already anticipated, the finite size effects on x"( w ) 
are quite weak. Both even and odd spectra show the 
presence of a prominent peak at a frequency of order of 
the gap A oc J to the first excited state. Note that in the 
absence of disorder, these peaks are Delta functions. The 
peak at wo ~ J for even systems results form excitations 
between states in the St o t=0 and St t=l sectors, while 
in the case of odd systems the transitions are between 
and 5(04=3/2 with an ujq ~ 3J/2 11]. However, 
we note that there is a fundamental difference between 
N even and odd in that the latter has significant low 
frequency spectral weight. In the rest of the paper, we 
shall discuss the behavior of this low frequency feature 
in the context of recent neutron scattering experiments 
on ZnCu 3 (OH) 6 Cl 2 0. 

The qualitative difference between the low energy spec- 
tra of even and odd systems can be readily understood. 
This spectra corresponds to transitions within the low- 
est energy sector that contains a macroscopic number of 
quasi-degenerate states (i.e., within energy 0(8 J) from 
the groundstate). This sector has St o t—0 for even sys- 
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FIG. 4: (Color Online) The low frequency peak in x"^) as 
a function of frequency uo in units of J for different temper- 
atures T= 0.001, 0.025 and 0.1 J. Taking J=170K, our data 
approximately correspond to the same experimental window 
reported in RefQ- The solid and dashed lines represent the 
power law fits proposed in that experimental study. 



terns and S't t=l/2 for odd systems. Since the spin sus- 
ceptibility induces transitions between states which dif- 
fer in one unit of the z component of the spin, these are 
forbidden between the singlets of even size systems, but 
are perfectly possible between the low-lying doublets of 
odd systems. Based on this argument, the low frequency 
feature in x"( w ) seen m °dd systems should have a char- 
acteristic width of order 5 J as is verified by Figs [3] and [4] 
Thus, the physics of the low frequency contribution to the 
dynamical susceptibility is also due to the same effective 
spin-i degree of freedom that was previously responsible 
for the low temperature contribution to x(T). 

We now focus on the temperature dependence of this 
low frequency peak. As one increases T from 10~ 3 J to 
1CP 1 J, which is roughly the same range as that studied in 
the experiments of Refi, we observe the following: first, 
the height and intensity of the peak diminishes dramati- 
cally and secondly, the peak shifts to higher frequencies. 
These behaviors are qualitatively similar to those seen in 
neutron experiments l 4!|. At smallest T there is a strong 
enhancement of x"( w ) towards lower frequencies. How- 
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FIG. 5: (Color Online) x"( w ) as a function of frequency ui. 
An external magnetic field is applied along the transverse x- 
axis (top) , the longitudinal 2-axis (bottom) and a tilted axis in 
the x — z plane (middle) . To match the experimental value we 
set h=0M7J, and T=0.001, 0.025 and 0.1 J (top, middle and 
bottom curve in each left panel). The left column panels show 
a detailed of the low energy region and the corresponding full 
spectra are shown in the right column (only the lowest T is 
plotted). 



ever, eventually the susceptibility vanishes as lo — > at 
the very low frequency end of the spectra. The curves 
for the three different system sizes shown in the figure 
demonstrate that even in this low frequency regime, the 
finite size effects are rather small. As argued before, the 
width of the peak is controlled by the disorder distri- 
bution 5 which is a phenomenological parameter in our 
model. Since J=170K, we adjust the disorder strength 
to the value S = 0.01 to fit the experimental neutron 
data. The frequency range shown in Figd] corresponds 
to the experimental range of the reported neutron scat- 
tering [4j and we clearly see that taking into account the 
experimental error bars, our results are indeed compati- 
ble with the observed data. For comparison, the power 
law fits suggested by the experiments are also shown in 
this figure. However, our results do not predict a diver- 
gence of f° r w ~* 0. Therefore, our study suggest 
that the putative divergence observed may be merely due 
to rather large error bars of the lowest frequency data, 
which is the region where the experiment suffers from 
largest uncertainty. Moreover, the lack of finite size ef- 
fects permits us to obtain a rather reliable estimate of the 
integrated intensity of the low frequency feature, and we 
find that it amounts to about 25% of the total spectral 
weight. This value is compatible with the reported 20% 
estimated from the neutron experiments. 

We finally consider the effect of a magnetic field on 
the dynamical susceptibility. In the presence of an ex- 
ternal field h = hz, x is the longitudinal susceptibility, 
whereas for fields aligned along the x— direction, x repre- 
sents the transverse susceptibility. In we plot the temper- 
ature evolution of x"( w ) f° r applied fields with different 
orientations. The strength of the magnetic field is set 



to h = 0.067J so as to match the experimental value of 
11.5Tesla used in Ref|4|. Our results show that h has a 
strong effect on x"( w ) at both high and low frequencies. 
For an applied field in the ^-direction, the longitudinal 
susceptibility x" (bottom panels in Fig(5j), remains es- 
sentially unchanged since a longitudinal field does not 
produce spin flips. For a transverse field, the computed 
x"(w) now corresponds to the transverse susceptibility 
and as expected, the high frequency feature at loq = 3/2 J 
shows the usual Zeeman splitting wo±/i (top right panel). 
This splitting corresponds to excitations from states with 
Sf ot = ±1/2 to states with Sf ot = ±3/2. Interestingly, 
we also see a shift of the low energy part of the spectrum 
towards a higher frequency lo ~ h. A single spin-^ would 
give a sharp peak at the Larmor frequency h. However, 
in the present case, in addition to a sharp peak, we also 
find a broad asymmetric contribution. This anomalous 
behavior is due to the quasi-degeneracy of the low-lying 
states and the many-body nature of the remnant spin- 
h degree of freedom present in odd systems. Interest- 
ingly, a similar broad anomalous feature was also seen in 
the reported neutron data For fields applied in the 
tilted (1,0,1) direction (central panels) the resulting re- 
sponse can be easily interpreted as the superposition of 
the transverse and longitudinal responses. Thus, the high 
energy feature is split into three peaks centered around 
loq and loq ± h, and a similar analysis is also valid for 
low frequencies. The response in a tilted field can be 
considered as a generic case, and is also more likely to 
be closer to the experiments which were performed on 
powder samples. 



We believe that the accord we find between our mean 
field toy model results and experiments may provide use- 
ful guidance, especially, in regard to the current debate 
on whether the low energy behavior seen in experiments 
is dominated by impurity effects. In fact, odd systems 
can be viewed as even systems with a spin defect (i.e., an 
impurity), or, alternatively, as a single (impurity) spin 
dressed by the coupling to an Stot=0 (even) spin liquid. 
Thus, our results for odd systems naturally lead to the 
interpretation of the anomalies observed in current ex- 
periments as arising from impurities/defects which may 
originate from the lack of precise stoichiometric control 
in currently available samples. Once a chemical handle 
is found for the synthesis of the compounds, a system- 
atic study of the dependence with Zn-Cu substitution 
may help disentangle the impurity and intrinsic contri- 
butions, and eventually validate our proposed scenario 
for this fascinating quantum magnet system. 



We thank C. Lhuillier for introducing us to this prob- 
lem. 
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